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Abstract 

We show that the Galactic latitude distribution of IceCube astrophysical neutrino events with energies above 100 TeV is 
inconsistent with the isotropic model of the astrophysical neutrino flux. Namely, the Galactic latitude distribution of the 
events shows an excess at low latitudes |6| < 10° and a deficit at high Galactic latitude |6| > 50°. We use Monte-Carlo 
simulations to show that the inconsistency of the isotropic signal model with the data is at > 3cr level, after the account 
of trial factors related to the choice of the low-energy threshold and Galactic latitude binning in our analysis. 

Keywords: multi-messenger astronomy; neutrino astronomy; Milky Way galaxy 


1. Introduction 

The discovery of astrophysical neutrino signal by Ice- 
Cube experiment has started a new field of neutrino as¬ 
tronomy The signal is a relatively soft 

(i.e. softer than dN^/dE oc E~^) powerlaw. Recent anal¬ 
ysis of the data @ excludes the l/E"^ type spectrum at 
more the 3.8 cr level. The signal is consistent with a 
powerlaw dN^/dE = A (A/100 TeV)”'"" with normalisa¬ 
tion A = 6.7/^fix 10“^® (GeV cm^ s sr)“^ and slope 

= 2.50 ±0.09 [3. 

It is not clear yet what kind of astronomical source(s) 
are detected by IceCube. The overall distribution of the 
three-year event sample is consistent with an isotropic dis¬ 
tribution [3, while a separate fit of the Northern and 
Southern hemisphere signals in the four-year signal shows 
a preference to a harder spectrum in the Northern hemi¬ 
sphere [3 which could potentially be due to the presence 
of a softer contribution of the flux from the inner Galaxy 
in the Southern hemisphere [3,13- Otherwise, the over¬ 
all approximate isotropy of the signal would point to its 
extragalactic origin. 

The unexpected properties of the observed signal chal¬ 
lenge pre-existing theoretical models. Softness of the spec¬ 
trum disfavours a range of models of extragalactic sources 
in which the neutrino flux is generated via interactions 
of high-energy protons with soft photon backgrounds. In 
such models charged pion production and decay which 
results in neutrino emission is characterised by a rela¬ 
tively high energy threshold. Neutrino energies in the 10- 
100 TeV range are typically much below the threshold and 
the neutrino spectrum in this energy ran ge i s expected to 
be much harder than observed 

Models which are favoured by the data are those in 
which the neutrino flux is produced in proton-proton (or, 
more generally, nuclei-nuclei) interactions. In this case. 


the energy threhshold for the pion production is as low as 
^ 100 MeV. The neutrino spectrum above this threshold 
approximately repeats the parent proton / nuclei spec¬ 


trum [1^. Measurement of T^ ~ 2.5 just tells that the 
parent proton / nuclei spectrum is a powerlaw with the 
slope Tp ~ Ti, ~ 2.5 Q- The parent protons / nuclei with 
such spectrum could be cosmic rays residing in the Milky 
Way galaxy (the Galactic component 

of neutrino flux). Otherwise, if the average spectrum of 
Galactic cosmic rays is much softer, the Galactic neutrino 
flux contribution to the IceCube signal is expected to be 
negligible [13 [13 and neutrinos have to originate from ex¬ 
tragalactic sources, like star forming galaxies [111 or 
radio galaxies / BL Lacs [13 113 (the extragalactic com¬ 
ponent of the flux). 

The extragalactic component of the neutrino flux is 
rather strongly constrained by the measurements of the 
isotropic y-ray background (IGRB) [2l|, Q. The IGRB 
spectrum is a powerlaw with the slope comparable to that 
of the neutrino spectrum, T.^ ~ 2.4, but with normalisa¬ 
tion which is approximately an order of magnitude lower 
than that of the neutrino spectrum 13 , 13 - At the same 
time, the fluxes of y-rays and neutrino prouduced in pp in¬ 
teractions are expected to have approximately equal both 
normalisations and slopes of the spectra 1J|- The IGRB 


constraint could be avoided if one assumes that the slope 
of the neutrino spectrum varies with energy, i.e. if the 
neutrino spectrum gets harder below ~ I TeV. The bro¬ 
ken powerlaw spectrum hardening at low energies is nat¬ 
urally expected in a situation when the low energy pro¬ 
tons are trapped in the source (or its host galaxy), while 
protons with energies higher than a certain threshold Er 
escape. This could be the case, e.g. for protons acceler¬ 
ated in the central engines and/or jets of radio loud active 
galactic nuclei [13 • The slope of the neutrino spectrum 
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changes around (0.01..0.!)_£,- [l^. Below this energy it 
follows original proton spectrum, which has power law in¬ 
dex 2.1 — 2.2 for usual Fermi acceleration models, while 
at higher energies diffusion of protons from the source 
through the turbulent magnetic field with Kolmogorov tur¬ 
bulence spectrum is expected to soften the spectrum to 
2.5 0. 

The possibility of non-negligible Galactic contribution 
is indicated by the consistency of the all-sky y-ray and neu¬ 
trino spectra, which follow the same powerlaw over some 
five decades in energy (from 10 GeV up to PeV) [3|. The 
y-ray all-sky spectrum is dominated by the Galactic con¬ 
tribution, so that it is natural to expect that the Galactic 
component is also present in the neutrino flux. The analy¬ 
sis of Ref. Q has searched for the correlation of the arrival 
directions of neutrinos with energies above 30 TeV in the 
three-year data set of IceGube. This analysis has found 
that the best correlation is at the level of ~ 2.5a pre-trial 
in the angle ±7.5° around the Galactic Plane and at the 
~ 2.2cr level (2.8% chance coincidence probability) after 
the trial factor is taken into account. 

Below we demonstrate that the neutrino four-year Ice- 
Cube signal in the energy band above 100 TeV iHi, 
which is free from the residual atmospheric neutrino and 
muon background Q, shows an evidence for the Galactic 
component. 

2. Anisotropy properties of neutrino signal 

The Galactic and extragalactic contributions to the 
neutrino flux could be distinguished based on the differ¬ 
ence in the expected distribution of the signal over the 
sky. The extragalactic flux should be isotropic, while the 
Galactic flux should show anisotropy toward the Galactic 
Plane, where most of the target material for the cosmic 
ray interactions is found. Low statistics of the neutrino 
signal and uncertainties in the modelling of the Galactic 
neutrino flux prevented a sensible analysis which would 
give definitive conclusions on the presence of the Galac¬ 
tic and extragalactic contributions in the first three years 
of IceGube data 0,113 ■ The overall distribution of neu¬ 
trino signal on the sky in the energy band above 30 TeV 
is consistent with an isotropic distribution ii , i.e. with 
the extragalactic signal. However, at 30 TeV the IceGube 
signal still has a significant contribution from the atmo¬ 
spheric neutrino and muon background which could dilute 
the weak anisotropy signal. 

A more clean anisotropy analysis could be performed in 
the energy band above 100 TeV, where the signal is almost 
backgorund-free. The updated results of 4-year IceGube 
exposure show 19 events in this energy band with only one 
background Fig. [T] shows the distribution of the 

detected E > 100 TeV neutrino events in Galactic latitude. 
One could notice two features in this distribution. First, 
the low Galactic latitude bin |6| < blow = 10° contains a 
large number of events (9 out of 19). Next, the bins at high 
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Figure 1: Galactic latitude profile of the E > 100 TeV IceGube 
neutrino signal. Dark grey solid histogram shows the expected profile 
of the isotropic neutrino signal. Dashed dark grey histogram shows 
the Galactic component profile. Thick light grey solid histogram 
shows the sum of the Galactic and extragalactic components. 

Galactic latitude (above bhigh = 50°) contain no events at 
all. 

To some extent, the lower number of counts in the bins 
at high Galactic latitude could be attributed to the smaller 
solid angle spanned by these bins. To verify if this would 
provide a satisfactory explanation of the deficit of neu¬ 
trino counts at high Galactic latitudes, we have performed 
Monte-Garlo (MC) simulations of the expected sky distri¬ 
bution of the neutrino signal. The MC simulation takes 
into account of the declination dependence of the IceGube 
effective area, derived from the information reported in 
the Ref. [3 . This declination dependence leads to a differ¬ 
ence in the effective exposure in different Galactic latitude 
bins. The MC simulation generates the number of events 
proportional to the exposure in each declination bin. The 
events are randomly distributed in the Right Accention. 
The MC events are then re-mapped in the Galactic coor¬ 
dinates. 

The Galactic latitude distribution of events expected 
in the isotropic flux model found from MC simulations is 
shown with the dark grey solid line histogram in Fig. [TJ 
If the isotropic flux is normalised on the total number of 
events, the isotropic model predicts 4.6 events at |6| > 50°. 
The probability to find no events in this latitude range is 
7'|b|>50° = 5 X 10 

The isotropic model is also inconsistent with the low 
Galactic latitude data, which shows an excess over the 
data. The tension between the model and the high Galac¬ 
tic latitude / low Galactic latitude data could be char- 
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acterised in a quantitative way using the MC simulations, 
which show that the probability to find simultaneously the 
0 counts at |6| > 50° and > 9 counts in the |6| < 10° bin 
is 10“^. This corresponds to a 4.4tT level inconsistency 
between the model and the data in terms of equivalent 
Gaussian statistics signal. 

The tension with the data could be readily removed via 
addition of the Galactic component of the flux. Following 
Ref. 0 , we model this Galactic contribution based on the 
E > 300 GeV y-ray data of Fermi/LAT. The Galactic com¬ 
ponent, convolved with the IceGube point spread function 
derived from the distributions of angular uncertainties of 
detected IceGube events, is shown by the dashed histogram 
in Fig. 1. One could see that the Galactic component ex¬ 
plains the excess in the first bin and it gives a negligible 
contribution at high Galactic latitude. The neutrino flux 
model which contains only Galactic component is also in¬ 
consistent with the data, but at a lower significance level 
(~ 2cr). A model which contains 50% contributions from 
the Galactic and extragalactic components provides a sat¬ 
isfactory fit to the data (Fig. 1). 

3. Account of the energy and angular binning trail 

factors 

The above conclusion on the 4.4tT level inconsistency of 
the isotropic model with the data adopts certain choice of 
the energy range and angular binning for the calculation 
of probabilities. To have a fair judgement of significance 
of this inconsistency, one needs to take into account the 
trial factor related to the choice of the Galactic latitude 
bins where the inconsistency of the model with the data is 
largest and to the restriction of the analysis to the events 
above a particular energy. 

There is no trial factor to be taken into account for the 
general choice of the pattern ’’high Galactic Latitude & 
low Galactic latitude” bins. Indeed, the highest / lowest 
Galactic latitude bins are those which contain the lowest 
/ highest Galactic signal. Any presence of the Galactic 
signal would reveal itself at the first place in these Galactic 
latitude bins as a deficit / excess of the data over the 
model. 

Typical width of the ’’low Galactic latitude” excess of 
the Galactic neutrino signal is determined by (a) the an¬ 
gular resolution of the IceGube telescope (which is about 
GpSF — 10° for the ” cascade” type events [Ei) and (b) 
by the angular width of the Galactic Plane itself. 

The width of the Galactic Plane depends on the di¬ 
rection on the sky. This ’’width” could be defined based 
on the y-ray observations of the diffuse Galactic emission 
reported in the Ref. Defining it as the angular dis¬ 

tance from the Galactic Plane at which the y-ray signal 
drops by a factor of 2, one could find that in almost all 
directions, the angular span of the Galactic Plane is within 
\b\ < 10°. This could be readily understood using the fol¬ 
lowing qualitative argument. Most of the target material 
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Figure 2: Top: dependence of the chance coincidence probability for 
the isotropic model to describe the E > 100 TeV data as a function of 
the cut defining the ’’low Galactic latitude” region (light blue 
curve on the left) and as a function of the cut defining the ’’high 
Galactic latitude” region (dark blue curve on the right). 


for cosmic ray interactions in the Galactic Plane is con¬ 
centrated in a narrow disk of the scale height H ^ 200 pc 
28, 29|. The density of the Disk drops beyond i? ~ 10 kpc 
28l . 1291 , This disk, observed from the position of the Solar 
system has different angular span in different directions. 
The estimate of the angular width of the Disk in the di¬ 
rection of the Outer Galaxy, 0 < H/{R — Rq) ~ 5° is 
consistent with the observed 0 ^ 10° extent of the pion 
decay y-ray signal in Galactic latitude in the outer Galaxy 
90° < I < 270° observed by Fermi/LAT [23|. The angu¬ 
lar span of the disk is smaller in the inner Galaxy region, 
shrinking to H/Rq ^ 1.5° in the direction of the Galactic 
Gentre [27|. 

Overall, the angular width of the Galactic component 
of the neutrino signal is expected to be somewhat larger 
than the average width of the Galactic Plane, 10°. This is 
illustrated in Fig. [1] which shows the Galactic latitude pro¬ 
file of the Galactic component calculated using the y-ray 
template convolved with the IceGube point spread func¬ 
tion for the cascade events. 

Light blue color curve on the left of Fig. [5] shows 
the dependence of the chance coincidence probability p 
on blow One could see that the probability is minimised 
to p cs 10“® at blow = 6.25°, which is comparable to the 
expected width of the Galactic component of the neutrino 
signal calculated based on the y-ray template. 

Dark blue color curve on the right of Fig. [2] shows 
the level of inconsistency of the isotropic model with the 
data as a function of the bhigh parameter. One could see 
that the chance coincidence probability reaches a broad 
minimum in the range 30° < bhigh < 60°. 

o take into account the trial factor introduced by choos¬ 
ing the values of blow, bhigh at which the inconsistency of 
the isotropic model with the data is maximal, we allow 
blow and bhigh to vary in the limits 0 < blow < b < 30° 
(as in the Ref. i) and 40 ^ bhigh ^ 90 . We esti 
mate the post-trial probability of the isotropic model to 
describe the data by counting the number of Monte-Garlo 
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Figure 3: Dependence of the chance coincidence probability for the 
isotropic model to describe the data as a function of the low energy 
threshold. 

simulated data sets in which the excess of event counts 
at |6| < blow and deficit of events at |6| > bhigh leads to 
the chance coincidence probability < 10“® . This results 
in the post-trial probability p = 2 x 10“^, or a 3.7cr level 
inconsistency of the model with the data. 

The inconsistency of the spatial distribution of the neu¬ 
trino signal with the isotropic model is expected to de¬ 
crease with the decrease of the energy threshold of the 
analysis. Indeed, at lower energies, the astrophysical sig¬ 
nal gets diluted by the increasingly important atmospheric 
muon and neutrino backgrounds. The background compo¬ 
nents obviously do not correlate with the Galactic Plane 
and their Galactic latitude distribution is similar to that of 
the isotropic astrophysical flux component. Thus, the de¬ 
viations from the isotropic distribution are more difficult 
to detect at low energies. 

Fig. |3] shows the dependence of the level of inconsis¬ 
tency of the isotropic flux model with the data as a func¬ 
tion of the low energy threshold. One could see that the 
chance coincidence probability reaches a minimum at the 
threshold E = 100 TeV. At lower energies the probabil¬ 
ity rapidly increases because of the increase of the back¬ 
ground. At higher energies the probability increases be¬ 
cause of the decrease of the signal statistics. 

We estimate the trial factor due to the particular choice 
of the energy threshold by extending the analysis down 
to 30 TeV ^he low energy threshold of the reported Ice- 
Gube data @,[1^). The number of generated Monte-Carlo 
events of different energies follows the energy distribution 
of the really observed events. The post-trial probability 
for the isotropic model to describe the data is calculated 
by counting the number of occurrences of the low Galactic 
latitude excesses and high Galactic latitude deficits lead¬ 
ing to the inconsistency level p < 10“^ in the simulated 
data sets after the adjustment of the energy threshold to 
the value which minimises the probability for the isotropic 
model to describe the data. This results in the post-trial 
probability p = 1.7 x 10“^, i.e. to the inconsistency of the 
isotropic model with the data at 3.1tT level. 


4. Comparison with the previous search of the Galac¬ 
tic component 

Search for the excess of neutrino signal in the direction 
of the Galactic Plane was previously reported in Ref. Q 
based on the analysis of the three-year dataset of IceGube. 
The difference between our approach and that of the Ref. 

0 is 

• the account of not only the low Galactic latitude 
excess, but also of the high Galactic latitude deficit, 
characteristic for the Galactic flux component and 

• dedicated analysis of the high-energy part of the neu¬ 
trino spectrum. 

As it is explained above, both modifications potentially 
lead to an increase of the sensitivity of the search via more 
precise account of the signal morphology and via suppres¬ 
sion of the unrelated background of atmospheric muons 
and neutrinos. 

To assess the improvement of sensitivity due to these 
modifications, we have repeated the analysis restricting the 
data set to the three-year IceCube data of Ref. Q . We find 
that the chance coincidence probability for the isotropic 
model to describe the three-year data above 100 TeV en¬ 
ergy threshold isp = 1.3x10“^ pre-trial andp = 7.8 x 10“^ 
post-trial. This should be compared to the p ~ l'^_/ 
p = 2.8% pre-trial / post-trial probability found in Ref. [Sj 
without an additional verification of the high Galactic lat¬ 
itude deficit, but using the likelihood analysis technique. 
The optimal width of the Galactic Plane is found to be 
blow = 6.25° (compare to the 7.5° found in the Ref. Q) 
and the optmal energy cut is at A = 200 TeV. 

5. Conclusions 

To summarise, we have shown that the Galactic lati¬ 
tude distribution of the neutrino events in the 4-year Ice- 
Cube data set at highest energies E > 100 TeV is inconsis¬ 
tent with the hypothesis of the isotropic neutrino flux at 
about 3(7 level. This inconsistency suggests the presence 
of a significant Galactic contribution to the signal. The 
statistics of the signal is still to low for the measurement 
of the relative importance of the Galactic and extragalactic 
flux components. A broad range of models with sizeable 
(the same order of magnitude) contributions of the two 
components are consistent with the present data. Tak¬ 
ing the 3 ct evidence level found in the 4-year data sample 
as an indication, one could estimate that the final verifi¬ 
cation of the presence of the Galactic component of the 
neutrino flux at A > 100 TeV would require a decade-long 
IceGube exposure. A much faster verification of the Galac¬ 
tic signal will be possible with the IceGube Generation 2 
detector [s^. Northern hemisphere neutrino telescopes, 
like km3net, should be able to detect the Galactic compo¬ 
nent of the signal at much lower energies due to the lower 
atmospheric background in the up-going muon neutrino 


4 





channel. The results of analysis of the data of ANTARES 
neutrino telescope show that already the sensitivity of this 
telescope is marginally reaching the expected flux level of 
the Galactic component (provided that it constitutes 100% 
of the IceCube signal) [3l|. This shows that a detailed 
study of the spectral and morphological properties of the 
Galactic signal will be possible with the km^ scale detec¬ 
tor, kmSnet [3^ . 
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